Chiral symmetry breaking, color superconductivity and quark matter phase diagram: 

a variational approach 
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We discuss in this note simultaneous existence of chiral symmetry breaking and color supercon- 
ductivity at finite temperature and density in a Nambu-Jona-Lasinio type model. The methodology 
involves an explicit construction of a variational ground state and minimisation of the thermody- 
namic potential. There exists solutions to the gap equations at finite densities with both quark 
antiquark as well as diquark condensates for the "ground" state. However, such a phase is thermo- 
dynamically unstable with the pressure being negative in this region. We also compute the equation 
of state, and obtain the structure of the phase diagram in the model. 
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I. INTRODUCTION 



The structure of vacuum in Quantum Chromodynamics (QCD) is one of the most interesting questions in strong 
interaction physics Jl| . The evidence for quark and gluon condensates in vacuum is a reflection of its complex nature 
H, whereas chiral symmetry breaking is an essential feature in the description of the low mass hadron properties. 
Due to the nonperturbative nature of QCD in this regime different effective models have been used to understand 
the nature of chiral symmetry breaking. These have been constructed, for the most part, in the framework of a 
Nambu-Jona- Lasinio (NJL) model with a four fermion interaction. There also have been attempts to generalise this 
£T) ' in the case of Coulomb gauge QCD with an effective propagator simulating the effects of confining potentials 

Studies at finite temperatures using imaginary time formulation of finite temperature field theory jll | have also been 
carried out. 

Recently there has been a lot of interest in strongly interacting matter at high densities. In particular, a color 
superconducting phase for it involving diquark condensates has been considered with a gap of about 100 MeV . The 
Q_i studies have been done with an effective four fermion interaction between quarks , direct instanton approach plj] 
or a perturbative QCD calculation at finite density There has also been a study of this phase in NJL model 

p5| . In the NJL model however ,the aspect of chiral symmetry breaking in presence of diquark condensates has not 
been considered. Such a question has been considered in Ref. [ p6| in an instanton induced four fermion interaction 
model within a mean field approximation. In this note, we use a different method to study the problem. We consider 
a variational approach with an explicit assumption for the ground state having both quark antiquark and diquark 
condensates. The actual calculations are carried out for the NJL model such that the minimisation of the free energy 
JH \ density determines which condensate will exist at what density. 

We organise this paper as follows. In section II we construct the ground state as a vacuum realignment with quark 
antiquark as well as diquark condensates. The finite temperature and density effects are included within the framework 
of thermofield dynamics JItJ through a thermal Bogoliubov transformation involving doubling of the Hilbert space. 
In section III we shall consider an effective model as is considered for chiral symmetry breaking in Coulomb gauge 
QCD f^-[|]. To solve the gap equation exactly we shall further take a four fermion point interaction limit of the same 
similar to NJL model. In section IV we solve the gap equations for the mass and the superconducting gap, determine 
the phase diagram and discuss the results. Finally, we give a summary in section V. 



II. AN ANSATZ FOR THE GROUND STATE 



As noted earlier we shall include here the effects of both chiral symmetry breaking as well as diquark pairing. For 
the consideration of chiral symmetry breaking, we shall take the perturbative vacuum state with chiral symmetry as 
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|0). In this basis we shall take the quarks as massless. We shall then assume a specific vacuum realignment which 
breaks chiral symmetry because of interaction. 

We have seen earlier that chiral symmetry breaking takes place with the formation of quark antiquark condensates 
in perturbative vacuum We consider the quark field operator expansion in momentum space given as 

"To||ilJT| 

[f/o(k)g?(k) + y (-k)g?(-k)] e lkx dk, (1) 
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(2tt) 3 / 2 
where pOpl 



The subscript I and the superscript indicate that the operators q® and q® are two component ones which annihilate 
or create quanta acting upon the perturbative or the chiral vacuum. 

We now consider vacuum destabilisation leading to chiral symmetry breaking (|JlTJ,|l^,|l?|| described by, 

\voc)=Uq\Q), (3) 

where 



U Q = cxp (J • k)h{k)q°j(k)dk - h. 



(4) 



In the above, we have suppresed the flavor and color indices on the two component quark and antiquark operators. 
Further h(k) is a real function of |k| which describes vacuum realignment for quarks of a given flavor. We consider 
here two flavors and take the condensate function h(k) to be the same for u and d quarks. Clearly, a nontrivial h(k) 
shall break chiral SU(2)l x SU(2)r symmetry to the custodial symmetry SU(2)y for the light quark doublet Jl9| , p2]| . 
In what follows we shall exploit this result. 

Having defined the state as in Eq.^) for chiral symmetry breaking, we shall next define the state involving diquarks. 
We note that as per BCS result such a state will be dynamically favored if there is an attractive interaction between 
the quarks p3] |. Such an interaction exists in QCD in the qq color antitriplet, Lorentz scalar and isospin singlet 
channel. In the flavor antisymmetric channel the interaction can be scalar, pseudoscalar or vector whereas in flavor 
symmetric channel only the axial vector channel is attractive. In the present work, we shall consider the ansatz state 
involving diquarks as 



where 

B 



d 
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|fi) = U d \vac) = exp(Bt - B d )\vac), (5) 



C(k)t/(k)^(-k)t^ e3afc + a*°(k)/i(k)er(-k)etfe3j dk. (6) 



In the above, i,j are flavor indices, a, b are the color indices and r(= ±1/2) is the spin index. As noted earlier we 
shall consider systems with two flavors and three colors. We have also introduced here two trial functions f(k) and 
/i(k) respectively for the diquark and diantiquark channel. As may be noted the state constructed in eq.(|3|) is spin 
singlet and is antisymmetric in color and flavor. The corresponding Bogoliubov transformation for the operators is 
given by 



C(k) \ = ( cos /( k ) -r F §SW e ^3 ac sm/(k) \/ g j»( k ) 

gff r (-k)t j ^ r rgi efcie3cQsin/(k) cos/(k) ) { gfi r (-k) 



(7) 



In a similar manner one can write down the Bogoliubov transformation for the antiquark operators corresponding to 
basis. 
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Finally, to include the effect of temperature and density we obtain the state at finite temperature and density 
|f2(/3, [/,)) by a thermal Bogoliubov transformation over the state |f2) using thermofield dynamics (TFD) as described 
in ref.s g^HH) • We have > 

\Sl{fi,n))=UpM (8) 
= e st (M- B (M, (9) 



where Up^ is 



with, 

B\(3,fx) = y^(k)t0_(k,/3, Ai )g;(k)t + g;(k)0 + (k,/3, M )^(k)dk. (10) 

In Eq.(|l0|) the ansatz functions 9±(k, (3, fi) will be related to quark and antiquark distributions and the underlined 
operators are the operators in the extended Hilbert space associated with thermal doubling in TFD method. Thus 
the ansatz functions at finite temperature and density given in Eq.(|^) invloves five functions - h(k), for the quark anti 
quark condensates, /(k) and /i(k) describing respectively the diquark and diantiquark condensates and 0±(k, (3, /i) to 
include the temperature and density effects. All these functions are to be obtained by minimising the thermodynamic 
potential. This will involve an assumption about the effective hamiltonian . We shall carry out this minimisation in 
the next section. 



III. ESTIMATION OF THE PRESSURE AND THE MODEL 



Since the low and medium energy behaviour of QCD is not well understood one has to make assumptions about the 
effective interaction between the quarks. To consider chiral symmetry breaking in the Coulomb gauge pairing model 
H-|| one had the effective action 

S eff = J d^ix^rd^ix) + \Jd 4 xJ dWi* ~ y)j;D^{x - y)J b u (y), (11) 

where, 

J^(x)=^(x) lfl (^ k (x). 

A sum over color indices a and b is implied, g is the strong coupling constant and, D^(x — y) is the full gluon 
propagator. This effective action may be thought of as resulting from integrating out the gluonic degrees of freedom 
from the full QCD Lagrangian keeping only the bilinear terms of the quark currents. While dealing with the gluon 
sector, to make things computable, one uses different approximations for the full propagator which is written down 
in Coulomb gauge [prpH^J ■ The most obvious one is to neglect transeverse gluons completely maintaining only the 
instantaneous Coulomb interactions. This was done in Ref. |i|,|5| for chiral symmetry breaking. Another possibilty 
is to neglect retardation effects for the transverse gluons and have e.g. D v (p) = _D u (p,po) - the so called Breit 
interaction. In such instantaneous approximations one can interpret the product of the coupling constant and the 
gluon propagator as the Fourier tranform of the effective quark antiquark potential. Different effective potentials such 
as confining, Richardson, and a screened potential for the transverse part have been considered for chiral symmetry 
breaking in Coulomb gauge pairing model ^-^| . 

To deal with the case of finite temperature and density one might take one loop resummed perturbative QCD 
results at finite temperature and densities |l4| , p7| ] for the propagator. However, we shall here assume a point like 
interaction approximation. The reason we shall choose this is its simplicity regarding its solvability. Further we can 
compare the results with those of NJL model calculations [jl5| . The delta function interaction produces short distance 
singularities and so to regulate the integrals we shall restrict the phase space inside the sphere |p| < A. Thus the 
effective Hamiltonian we shall be considering is given by 

a 2 

n = t/jH-ia- v)v + y j Z j>m - ( 12 ) 

We next write down the expectation values of various operators in the thermal vacuum given in Eq.(^). These 
expressions would be used to calculate thermal expectation value of the Hamiltonian to compute the thermodynamic 
potential. We define 
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and, 



where, 



and, 



(fi(/3, M )|^(k)^(k') t l^(/3, M)> - A +Q(3 (k, 0, M )<5(k - k'), 
(Sltf, M )|^(k)t^(k')|0(/3, m)) = A_ a/3 (k, /?, M )<5(k - k'), 



A+(k, (3, fi) = -(l + Fi (k) - F(k) + ( 7 ° sin 2h(k) + a ■ k cos 2/i(k)) (l - F(k) - F a (k))j , 



A_ (k, /3, = r (l + F ( k ) - F i ( k ) - (7° sin 2h(k) + a-k cos 2/i(k)) (l - F(k) - Ft (k)) 



(13) 
(14) 

(15) 
(16) 



Here, the effect of diquark condensate and the temperature and/or density dependence is encoded in the function 
F(k) and Fi(k) given as 



and, 



F(k) = sin 2 /(k) +cos2/(k) sin 2 6L(k, (3, fi), 



Ft (k) = sin 2 h (k) + cos 2ft (k) sin 2 9+ (k, (3, fi) . 



(17) 



(18) 



Clearly, at zero temperature and zero density the functions F and F± vanish and the projection operators reduce to 
the forms considered earlier ||lOf , 



A±(fc,/3) = i 



1 ± (7 sw2h(k) + a ■ k cos2h{k)) 



Further, at zero density but finite temperature 6>_ = 9 + and projection operators reduce to p5|, 



A±(£,/?) = ~ 



1 ± cos 26»( 7 ° sm2h{k) + a-k cos2h(k)) 



(19) 



(20) 



We also have 



W,fi)\^(xW^{P,fi)) = J e^V +7a (k,P,fi)dk, 

J e tk x V^ ai (k,(3,fi)dk, 



W,^\xW^W,n)) 



e ij e 3ab 

~(2nf 



where, 



V+ (k, f3, = [S(k) + (7° sin 2h(k) - a • k cos 2/i(k)J A(k) 



75C, 



and, 



V- (k, (3,fi) = (S(k) + (7 sin 2h{k) - a ■ k cos 2/i(k)) A(k) 



(21) 



(22) 



(23) 



Here, C — z 7 2 7 ° is the charge conjugation matrix (we use the notation of Bjorken and Drell) and the functions S(k) 
and A(k) are given as, 



S(k) = sin2/(k) cos 26L(k, /3, fi) + sin2/i(k) cos26» + (k, (3, fi), 



and, 



A(k) = sm2f(k) cos 29 -(k, (3, fi) ~ sin2/i (k) cos 29+ (k, (3, fi) , 
Using Eq. (|l4|) we have for the kinetic energy part 



(24a) 
(24b) 
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T=(Q,^fj,\^(-ia- V)V|W = y<flccos2fc(k)(l-F-Fi) ) 



(25) 



where, the degenracy factor 7 = 12 corresponding to spin(2), color(3) and flavor(2) degrees of freedom. T in Eq.([25|) 
above, however, includes the contribution of the perturbative zero point energy corresponding to h(k) = = /(kj = 
fi (k) . Subtracting this out we have the contribution from the kinetic energy 



T = T - (T)|/ l= / =/l=9±=0 
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(2tt): 



/ |k| [2sin 2 /i(k) +cos2/i(k)(F + F 1 ) 



(26) 



Similarly the contribution from the interaction term in Eq.(|l2]) after subtracting out the zero point pertuturbative 
energy turns out to be 



V = (fi, /?, M |y J° J"°|W = + ^2 
where, the contribution Vi arises from using Eq.s (|l3|), (n4h and is given as 



Vi = \N f 



(N 2 C - 1) 



2(1? - 2/f) = Sr? 2 ^ 2 - 2/f) 



with, 



and 



72 = (2^)3 / dk ( 1 -- F "-- F i) sin2 ' l (k)- 



(27T) 

The term V2 arises from using Eq.(22)and ( p3| ) and we have 

K 2 = ^(-2/| + /I), 



where 



and 



(2tt) 3 



dkS(k) 



(2tt) 3 



cflo4(k) sin2/i(k). 



Then the free energy density is 



T = e- nN = T + V - nN 
where, [i is the chemical potential and N is quark number density. Further, 



Finally, for the entropy density we have |l 

7 



{2-Kf 



dk (sin 9 -In sin + cos 9 -In cos + sin 9 + lnsin 9+ + cos cos 9^ 



(27) 



(28) 



(29) 



(30) 



(31) 



(32) 



(33) 



(34) 



(35) 



(36) 



and the thermodynamic potential which is negative of pressure is given by [|26| 



fi = -V = e - fiN - is 



(37) 
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Now if we minimise the thermodynamic potential with respect to /i(k),we get 



WMk) = * ^ (38) 
where, M = 8g 2 l2/3. Substituting this back in Eq.(|30|) we have the mass gap equation 

3 3 (2tt) 3 7 Vk 2 + M 2 V U 1 ' 

Clearly, the above includes the effect of diquark condensates as well as temperature and density through the functions 
F and F\ given in Eq.s flL7| ) and ( p"8| ) respectively. The zero temperature and density limit is given by setting 
F = = Fi which has the same structure as in NJL model fi ^RM . 

Next, minimising the thermodynamic potential Q with respect to / and /i yields 

, ofr .s 2g 2 (2/3-/4 sin 2h) 

tan2f(k) = — ? - x (40) 

9 Vk 2 + M 2 - (ji - ifh 

and, 

to n2/ l( k) = 2 4 Qh + h^h) r . (41) 



9 Vk 2 + M 2 + (/i- 



We might further simplify equations ( |40| ) and ( M]J ) by noting from Eq.s (|32|), (|33|), (24b) and (24a) that the integral 
Ii is small compared to ^3, as the integrand in the former is a difference of two postive quantities whereas the latter 
is a sum of these two quantities. In addition, the integrand in J4 is suppresed by the quark antiquark condensate 
function sin2/i(k) which is not the case in I3. Finally, this apart, in Eq.s (ffOj) and ( f4l|) the numerical coefficient in 
the second term of the numerator is small compared to that in the first term. Thus, we can approximate 

tan2/(k) = (42) 



ton2A(k) = - A-, (43) 



where, the superconducting gap A = (4<? 2 /9)i3, E = Vk 2 + M 2 and v = fj, — 4g 2 /3ii is the chemical potential in 
presence of interaction |fT^| . From the definition of superconducting gap and Eq.(g2|) we have the superconducting 
gap equation given by 

A = ^—h = [ dk I A = cos 26L (k, p,p)+ A = cos 29+ (k, 0, /x) ) (44) 

9 9(27r) 3 y \^/A 2 + (E-v) 2 V ; ^/A 2 + (E + v) 2 V 'J V ' 

Finally, minimisation of the thermodynamic potential with respect to the thermal functions Q± (k, /?, /i) gives 

sin 2 6L = (45) 

cxp(/3w_) + l v ' 

and 

sin 2 tf+ = 1 (46) 

exp(/^o; + ) + 1 



where, w± = y A 2 + £+ and £± — (E±v). Therefore, with all the ansatz functions determined from the minimisation 

of the thermodynamic potential fl(/3,fi), the mass gap equation (|3^) and superconducting gap equations are given 
respectively as 



rfk —tanh{- ) + — tanh(—— ) = 1 (47) 



3 (2tt) 3 J Vk 2 + M 2 \u- y 2 ' lu+ 2 
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' ^f-fei.-iML, (48) 



9 (2tt) 3 J \ lu_ uj + 

Equation (^) is the generalisation of the mass gap equation considered at zero temperature and density of Ref. 
|7|,0[ to include the effect of finite temperature and density [ [l5[ alonqwith the effect of diquark condensates. Similarly, 
Eq. fl48|) is the relativistic generalisation of the BCS gap equation |L5| in the presence of a dynamically generated mass 
gap through a quark antiquark condensate structure for the vacuum. These arc the main new features of the present 
work. 

IV. SOLUTION OF THE GAP EQUATIONS AND RESULTS 

Before trying to solve the coupled gap equations ( |47| ) and ( |4S| ) we discuss the solutions in the following two different 
limiting cases, so that we make connection with earlier studies JT^JI^]. 

• Case-I :Only quark antiquark condensate i.e. h(k) ^ 0; /(k) = = /i(k) 
In this limit A = and the gap equation (|47]) reduces to 

M= 8 f / . M =(l-n + -n-)dk (49) 

3 (2tt) 3 J Vk 2 + M 2 V + ' K ' 

where, 

1 

n± 



exp(/3(£(k) ±v)) + l 



This is the same equation as obtained in Ref. g,|28|,29j except for numerical factors before the integrand. This 
is due to the difference in the Lorentz structure of the four fermion interaction term. In Ref. fjv the interaction 
term is (4>ip) 2 + {{'4>lb T i } ) 2 that can originate from a Fierz transformation of the J° J Ma interaction as considered 
here. In Ref. j?9) only Jq J 0a term was included. We might however mention here that keeping Jo Jo term only 
will make the pressure negative at high densities even when the condensate vanishes. At zero temperature the 
mass gap is given by the solution of the equation 



,2 



k 2 dk 



Sir 2 J kf V^ 2 + M 2 



= 1 (50) 



where the fermi mometum fc/is defined through the (interacting) chemical potential v as kj + M = v . We 
also note that this gap equation is derived here by minimising the thermodynamic potential over a variational 
ansatz and not through a mean field approximation p§] or a Hartree Fock calculation Q. The critical fermi 
momentum is given by the density where M = and is given at zero temperature by, 



fc ^ A < 1 -w ,,/J (51) 



and the corresponding quark number density is given by no = 7fc^ 3 /67r 2 . The variation of mass gap as a 
function of fermi momentum is shown in Fig. (|l|). We have chosen the value of the coupling g 2 — 55.98GeV~ 2 
and A = 0.67GeV as typical values. In fact, they have been chosen so that the diquark condensate function 
vanishes at the same critical temperature as in Ref. 
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FIG. 1. Mass gap vs fermi mornetum. Critical fermi mometum turns out to be 0.39 GeV corresponding to a quark number- 
density of 1.7 per fm 3 . 



With this choice of parameters, the dynamically generated quark mass at zero temperature and density becomes 
about 490 GeV which is rather high compared to the standard value of around 300 MeV. With a lower value 
of g 2 A 2 the same could be made smaller. The critical density in this model turns out to be about 1.7/ fm 3 
corresponding to a value of the fermi momentum k'j of OAGeV. In Fig. 2 we have plotted the pressure as a 
function of fermi momentum. Here we have added the bag constant eo which is the energy density at zero quark 
number density at zero temperature |||. This turns out to be about — (0.2GeV) 4 . As observed earlier |Q the 
pressure has a cusp like structure and becomes negative at finite density. The portion of the curve that goes 
down with fermi momentum corresponds to a nontrivial mass gap solution and the portion that increases with 
density correponds to zero mass solution of the mass gap equation (f47j). As in Ref |jL5|| , and unlike Ref. JlJ] the 
effects of interactions do not disappear beyond the chiral symmetry restoration phase. This is due to the fact 
that there is no running coupling involved here and in fact, the pressure grows as n 2 or k® at high densities as 
may be expected from the interaction term V\ involving I 2 = N 2 /j 2 given in Eq.(28). However, one has to keep 
in mind that to be consistent with the philosophy of high momentum cutoff A, the fermi momentum should be 
less than the cut off momentum. The negative pressure at intermediate densities can be understood in terms 
of mechanical instability and can have the interpretation that uniform nonzero density quark matter will break 
up into droplets of finite density in which chiral symmetry is restored surrounded by empty space with zero 
pressure and density. It is tempting to identify these droplets with nucleons within which the density is nonzero 
and (qq) =0 — a fact reminiscent of bag models. Nothing within the model however implies that the droplets 
have quark number three |B3. 
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FIG. 2. Pressure as a function of fermi mometum. The decreasing part of the curve corresponds iom/0 solution and the 
rising part of the curve corresponds to m — solution of the mass gap equation (UjO 



Case- II:No chiral condensate but only diquark condensate i.e. h(k) = /(k) ^ 0, /i(k) ^ 



In this case, the quark mass arising from chiral condensate is zero and the superconducting gap equation (|4 
reduces to 



2.9 



2 r A 



k A dk 



tanh(— ) H tanh{ ^ + ) 

uj- 2 o;+ 2 



(52) 



where, as before u± = y A 2 
equation in superconductivity 



£4 but, with £± = |k| ± v. This is the relativistic generalisation of BCS gap 

3) . It has the same structure as in Ref in the limit the form factor is replaced 
by a sharp cutoff as in Ref. |15] and apart from the numerical factor preceeding the integrand. Further, this can 
be rearranged to interpret separately contributions arising from particles, antiparticles and holes |12|. We have 
chosen the values of the coupling and the cutoff as in Ref. |l5|] so as to get the same critical temperature as in 
Ref. p2[ . The gap equation at zero temperature is given as 



1 = 



9tt 2 



,A 


1 


1 " 


/ k 2 dk 






'0 




UJ + _ 



(53) 



The superconducting gap A is plotted in Fig ||. While obtaining the solution of the gap equation (^) we also 
insist that a nontrivial solution of the gap equation is acceptable only if the corresponding free energy is smaller 
than that with no gap. The gap increases with fermi momentum, has a maximum of about 90 MeV around 
fermi momentum 550 MeV, beyond which the effect of cutoff is felt and it vanishes at around kf = 600 MeV. 
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FIG. 3. Superconducting gap vs fermi mometum. 

The resulting equation of state i.e. pressure as a function of fermi momentum is plotted in Fig. 4. There is no 
unstable phase as discussed for chiral condensate case with negative pressure. Because of the interaction term, 
the equation of state does not go over to the free massless equation of state beyond the critical density. We have 
plotted the pressure in Fig. 4 for fermi momentum upto 600 MeV noting that the cutoff here is about 650 MeV. 
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FIG. 4. Pressure vs fermi mometum. 



Case- III: h(k) ^ 0;/(k) ^ Oand /i(k) ^ 



Here we solve the coupled gap equations (47) and (|48|) which at zero temperature but finite density reduce to 
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2g 2 


f d 4 






37T 2 J 






UJ+J 



= 1 



(54) 



9tt 2 



dkk 2 [ — + — )=1 
- ^+ 



(55) 



where, w± = w A 2 + £± = E±v and £7 = Vk 2 + M 2 . To solve the above equations numerically, we take v 

as an input and obtain the values of mass gap M and the superconducting gap A so that both the equations are 
satisfied simultaneously. The chemical potential v is related to the fermi momentum as kt = yjv 2 — M 2 . The 
resulting mass gap and the superconducting gap are plotted in Fig. 5 and Fig 6 respectively. For the sake of 
comparison we have also plotted the mass gap without diquark condensate (case-I) in Fig 5, and, superconducting 
gap without chiral gap (case-II) in Fig. 6. As may be noted, the mass gap does not change very much through 
the inclusion of diquark condensates. The superconducting gap however starts at a lower threshold and varies 
slowly compared to the case of no mass gap. Both the curves however merge at the chiral restoration point. 
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FIG. 5. Mass gap M as a function of fermi mometum. The dashed line correspond to no 
corresponds to both diquark and quark antiquark condensates. 



condensates. The solid line 
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FIG. 6. Superconducting gap A as a function of fermi mometum . The solid line corresponds to both diquark and quark 
antiquark condensates. The dashed line corresponds to only diquark condensates. 



We have also plotted the equation of state , pressure as a function of fermi momentum in Fig.^. The equation 
of state does not change very much compared to the case with only quark antiquark condensates. 
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FIG. 7. Pressure as a function of fermi momentum . 



This should be expected as the effect of diquark condensate is small in the region where chiral condensate 
is nonvanishing. As may be evident from Fig.^j) there is a region from kf ~ 0.1 fm till kf ~ 0.4 fm where 
simultaneous existence of both types of condensates is possible. However, in this region the pressure becomes 
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negative. As in case I, this will correspond to mechanical instability. The finite density matter will break into 
droplets of finite density in which chiral symmetry is restored and the matter is in a superconducting phase 
surrounded by empty space with zero pressure and density. Below kf = 0.1 fm the phase with only quark 
antiquark condensates and for kf > 0.4 fm the phase with diquark condensates is thcrmodynamically feasible. 
However, one should not extrapolate to fermi momenta higher than the cutoff momentum of the model to be 
consistent with the philosophy of high momentum cutoff. 



A. The phase diagram 

To discuss the phase diagram we have solved the two gap equations Eq. ( JT^ ) and Eq. ( [48| ) at finite temperature and 
density to calculate different thermodynamic quantities. The result of such a calculation for the mass gap is shown 
in Fig.(||). Similarly the superconducting gap at different temperatures is shown in Fig.(||). 




FIG. 8. Mass gap at different temperatures as a function of fermi momentum. 

The superconducting gap decreases with temperature and vanishes at around 55 MeV. Similarly the mass gap also 
decreases with temperature and vanishes at about 270 MeV at zero density. The superconducting transition is a 
second order phase transition. The chiral phase transition is first order at zero density but is second order at high 
temperature. The pressure at different temperatures is shown in Fig. 10 as a function of (n/no) 1 ^ 3 , where n is number 
density and no = (2/tt 2 )v^ is the critical number density at zero temperature. The cusp in the pressure density 
curve, indicating a first order phase transition, vanishes at about 84 MeV. At temperaures beyond that the pressure 
increases monotonically with density. 
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FIG. 9. Superconducting gap at different temperatures as a function of fermi momentum. 
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FIG. 10. Pressure as a function of density at different temperatures. The tricritical point turns out to be about 84 MeV. 



To discuss the critical line for the chiral transition in the number density and temperature plane we need to solve 
the gap equation ( [l7| ) with mass gap equated to zero i.e. the critical line satisfies the equation 

dkkil T^TTT 1 -T7T 1 ~ T ) = 1 (56) 



4g 2 



/0 v exp(f3(k - v)) + 1 exp(/3(fc + v)) + 1, 

As the distribution functions are rapidly decreasing functions of momentum we may approximate the upper limit 
A — > oo. In that case Eq.(p39) reduces to 
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2,2 2 
V + T =V C 



(57) 



where, r = TtT/\fi and v c = k c f = A(l — 3ir 2 / (2g 2 A 2 ), the critical chemical potential or fermi momentum at zero 



temperature calculated earlierin Eq.(51). It may be useful to write down this relation in terms of number densities 



rather than the chemical potential using the relation 



fA dkk 2 ( — 1 .. * \k\ v (i? +1 ?1*) (58) 

Vexp(/3(fc - v)) + 1 exp(/3(fc + v)) + 1 / tt 2 v ; v ' 



n 



where, we again take infinity as the upper limit in the momentum integration as was done in Eq. (|5q) . This leads to 
the equation for the critical line in the number density plane as 



T 



1/2 



where, n = (2/tt 2 )v^ is the critical density at zero temperature. This approximate solution is shown by the dotted 
line in FigpT|. The numerical solution for the critical line with the finite cutoff for the momentum is shown by the 
solid curve in the same figure. As the chiral phase transition is a first order phase transition at zero temperature 
and is second order at high temperature there will be a tricritical point for the chiral phase transition below which 
there will be a mixed phase. This mixed phase here will correspond to the phase with droplets of finite density quark 
matter in the superconducting but chiral symmetric phase surrounded by empty space with n = 0, P = 0, A = and 
chiral symmetry broken phase. This can be obtained by a Maxwell construction to exclude the part of the phase 
diagram where pressure decreases with density. This is shown as the dashed curve in Figjll]. The tricritical point in 
this model turns out to be 84 MeV 

We have also shown in Fig.|ll| the critical line for the color superconducting phase by the dot-dashed line. The 
decreasing part of this line is due to the fact that we have taken a cut off in the momentum and at high enough 
density the superconducting gap vanishes as shown in Fig. [)| 




FIG. 11. Phase diagram in density temperature plane 

A detailed quantitative comparison of our results with those in Ref. |l6f| is not very meaningful since our model is 
a simple, schematic one. However, it is very reassuring that the structure of the phase diagram is very similar to the 
one obtained by Berges and Rajagopal fL6[ . 
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V. SUMMARY 



We have analysed in a current current point interaction model the structure of vacuum with quark antiquark as well 
as diquark pairs. The methodology used is a variational one with an explicit construct of the trial state. The present 
work is not based on a mean field approximation p8[ . Because of the point interaction structure we could solve for 
the gap functions explicitly. If we had taken the interaction term with a potential the gap equation would become an 
integral equation. It will be interesting to include a realistic effective potential and solve for the gap equation. There 
appears to be a region depending upon the coupling where both chiral condensates as well as diquark condensates 
are thermodynamically feasible. Presence of diquark condensate does not modify the dynamical mass of the quarks. 
The dynamical mass however affects the threshold for superconducting gap. The equation of state does not differ 
very much with inclusion of diquark condensates. We also obtain the complete phase diagram of the system in overall 
agreement with that of Ref. jl6| . Finally, our results suggests that our simple model contains the essential physics of 
the system. 
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